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ABSTRACT 

We report spectra from 3000-5900 A for comet 96P/Machholz, obtained on 
2007 Apr 27 UT with the 3-m Shane telescope at Lick Observatory. The spectra 
are extremely carbon poor, and show a prominent NH2 series, but no CN emis- 
sion. NH, NH2, and C2 gas production rates are (8.36 ± 2.18) xlO^^, (29.88 ± 
3.66) X 10^^, and (4.52 ± 0.61) x 10^^ molecules sec~^, respectively, as determined 
from Haser model fits to the data. Upper limits to the gas production rates for 
CN and C3 are 7.5x10^^ and 2.0x10^^, respectively. Though 96P is depleted in 
C2 and C3 relative to NH, it is even more depleted in CN than other so-called 
"carbon-chain depleted" comets. 



Subject headings: comets: general - 
Cloud — solar system: formation 



comets: individual (96P/Machholz) 
- (ISM:) cosmic rays 



Oort 



1. INTRODUCTION 

In recent decades a rapidly growing number of periodic comets have been discovered, 
including new members of the low-inclination short-period Jupiter family, and the high- 
inclination medium-perio d Halley family. A sub-cla s s of Jupiter farnily comets has been 



shown by many authors (INewburn fc Spinrad 1 1 19841 : lA'Hearn et al. I Il995l : iFink fc Hicks 



19961 ) to be depleted in the "carbon-chain" molecules C2 and C3 relative to CN. A few comets 
with more abnormal produ ction rate r a tios h ave been reported, including notably C/1998 
Yanaka (1988r; Fink 1992). ISchleicher I (120071 ) (hereafter S07) recently reported unusual gas 
production rate ratios for 96P, determined from narrowband photometry performed on 2007 
May 12 and 24 when the comet was at r^ = 1.07 and 1.30 AU, which we now confirm. 



Comet 96P/Machholz was discovered in May 1986 (IMachholz. Morris. &: Hale Ill986l ). 
Surprisingly short-period (5.2 years) orbital elements were calculated from th e initial magni- 
tudes and positions reported in the lAU circulars (Tabled]). ISekanina I (Il990al ) noted a major 
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outburst several weeks after perihelion and suggested that the comet remained undiscovered 
in spite of its short period because it is inactive for m uch of the time. Its small perihelion (0.12 



AU) leads to a high erosion rate (ISekanina Ill990bl ). so it may have a mostly dormant iceless 



about 3 AU from the Sun in 1997, and observed by SOHO on three occasions ( 


Meech et al. 


1997; 


Lisse, Fernandez. & Biesecker 


2002: 


Grvnko, Jockers. & Schwenn 


2004j). 



96P is d y nami cally interesting. It is the shortest-period Halley family comet known. 
Green et al. I (jl990l ) noted that 96P (then known as comet 1986VIII) travels closer to the 
Sun than all other known comets with orbital periods less than 150 years. It is in a 9/4 
resonance with Jupiter, and their orbit evolution calculations show that its perihelion is 
steadily decreasing to a minimum of 0.03 AU in about the year 2450. 96P has also been 



suggested as the parent body for several meteor streams (ISekanina Ill990al : iMcIntosh Ill990 



Babadzhanov fc Obrubov Ml993r). All other references to 96P discuss i ts meteor stream 



connection (jjenniskens et al. 1 119971 : IWilliams fc CoUander- Brown I Il998l ). and the possi 
bility that some sunskirting comets are fragments of it that were thrust onto their cur 



rent parabolic orbits after encounters with Jupiter (lOhtsuka. Nakano. &: Yoshikawa 1 12003 
Sekanina fc Chodas II2005I ). 



2. OBSERVATIONS 

Comet 96P/Machholz was observed post-perihelion on 2007 April 27 UT with the Kast 
Spectrometer at the Lick Observatory 3-m Shane telescope. The blue side was used with a 
1.5x128 arcsec slit and a 452/3306 grism giving a dispersion of 2.49 A/pixel. A Reticon CCD 
recorded the 2-D spectra (1199 pixels along the dispersion direction and 164 pixels along 
the slit at a scale of 0.78 arcsec/pixel). The comet's orbital parameters and observational 
circumstances are given in Table [H 

The data reduction used standard IRAF0 techniques and custom IDL routines. In IRAF, 
the spectra were trimmed to encompass the length of the slit, fiat-field corrected, median 
filtered to remove cosmic rays, and wavelength calibrated using HeHgCd lamp reference 
spectra. 

Because 96P rose shortly before morning astronomical twilight, there was time for one 
dark-sky comet frame to be taken starting at 11:38 UT, but not enough time to obtain a 
separate sky spectrum for background subtraction. The second comet spectrum was taken 



http: / /iraf. noao.edu/ 
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mostly between astronomical and nautical twilight and contained a large scattered sky con- 
tinuum component. Fortunately, both 2-D spectra contain sky information at the ends of 
the slit, with only a small amount of flux contribution from the NH band at 3360 A. In IDL, 
we constructed a virtual 2-D sky frame for each comet exposure by fitting a line under the 
NH flux at both the top and bottom of the slit to remove it, and linearly interpolating the 
sky values from the top and bottom rows of the CCD frame across the entire slit length. 
Back in IRAF, each virtual sky frame was subtracted from its comet frame. 

The final flux calibration was done with a custom IDL code that applies an extinction 
curve measured for that night, and utilizes standard fluxes blueward of 3200A not available 
with IRAF. The spectrograph range was set this blue to observe the OH 3085A band in other 
target comets. The extinction curve was constructed from observations of the standard star 
BD-l-33 2642 over the airmass range 2.52-1.12, using a wide (9 arcsec) slit to capture all of 
the starlight. The extinction versus wavelength relation that was derived is shown in Fig. [1] 

The IDL fiux-calibration code does several things: integrates a spectrum of the standard 
star (taken with the same 1.5 arcsec slit as the comet) along the slit, corrects the integrated 
standard star spectrum to zero airmass, and creates a fiux-per-count "sensitivity function" 
by dividing the standard fluxes by the integrated standard star counts for each pixel along the 
dispersion direction. The standard fluxes were taken from two sources: the magnitudes per 
unit frequency distr ibuted with IRAF redward of 3200 A, and the flux per unit wavelength 



from iBohlin I ( 19861 ) blueward of 3200A. This "sensitivity function" takes into account total 
system throughput, including wavelength- dependent light losses from using a narrower slit. 
The last IDL task is to calculate the appropriate extinction correction for each comet frame, 
and fold the comet extinction correction into the sensitivity function to create a set of final 
flux-calibration masks. Finally, back in IRAF, each 2-D comet spectrum is multiplied by its 
flux-calibration mask. 

The average of our two flux-calibrated 1200-sec spectra of comet 96P, integrated along 
the slit, is shown in Fig. [2l The NH emission band at 3360 A is clearly visible in Fig. [2^. 
Various NH2 emission bands are identified in Fig. [2)d. Most of the identified lines are NH2 
and NH, presumed to be photodissociation daughter species of ammonia. Notable is the 
faint C2 emission, and the questionable CN and C3 emission. Some flux was apparent in the 
region of the C3 band for one spectrum, but the other spectrum only showed noise in this 
location. 



-4- 



MOLECULE PRODUCTION RATES 



Molecule production rates were calculated from a iHaser I (119571) model. For a detaile d 
explanation of how this model w as applied to KAST data, please seelTurner fc Smith I (Il996[). 
Fluore scence efficiencies ar e fromlKim. A'Hearn. fc Cochran I (119891 ) for NH. lNewburn fc Spinrad 



(Il989h for C2 and C. andEatiim I tdsj) for CN. For NH^, we used the blue band fluores- 



cence efficiencies of iTegler fc Wvckoff I ( 



1989), scaled to the newer red band efficiencies of 



Kawakita fc Watanabe I (120021 ) . Following lA'Hearn et al. I (119951 ). the gas outflow velocity 
scales with heliocentric distance as 1.0 x r/i (AU)~°'^ km s~^. Sets of scale lengths were de- 
termined for the NH, NH2, and C2 bands in each spectrum by experimenting with a range of 
possible values. Where the best-flt scale lengths for NH2 bands blueward of (0,12,0) noisily 
exceeded 100x10^ km, the mean of the other c leaner scale lengths was used instead. Scale 
lengths for the CN and C3 bands followed from [Turner fc Smith I (119961 ). 



The scale lengths (applicable at the comet's heliocentric distance Vh = 0.754 AU), and 
the resulting gas production rates [Q) with 1-a error bars, are given in Table [21 "NA" in 
Table [2] refers to an undeflned production rate where the net flux in a non-detected emission 
band was less than the adjacent continuum. S07 report log Q(OH) = 27.33, while we did 
not detect OH. Our production rate for NH is a factor of 3 higher than S07, consistent with 
our smaller Vh than S07, but our production rate for C2 is 90% of S07's. Our upper limits 
for CN and C3 are consistent with S07. 96P may be sputtering somewhat unevenly as its 
global gas production turns off. 



4. DISCUSSION 

The only other comet in the literature noted to be this deflcient in CN and C2 is Yanaka 
^1988r or C/1998 Yl; Kosai et al. 1988). Yanaka was observed post-perihelion in 198 9 with 



the f/1.2 s 



1991 



Dectrograph at the University of Arizona 1.54-m Mount Bigelow telescope ( 



I992I ) . The observational parameters are listed in Table [3l where not given in 



Fink 



Fink 



(Il992h . they were looked up at the Minor Planet & Comet Ephemeris Service! provided by 



the lAU. 

Comet Yanaka displays a series of NH2 emission bands from 5300-8500A. Fink's spec- 
trum of it and our spectrum of 96P/Machholz overlap from 5300-5800A, in the region of the 
NH2 (0,10,0) and (0,11,0) bands. Noticeably absent in the Yanaka spectrum are the CN 1-0, 
2-0, 2-1, and 3-1 bands, and the prominent C2 Av = —1 band. 



2http://www.cfa.harvard.edu/iau/MPEph/MPEph.htnil 
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Greenberg. Singh. &: de Almeida I (Il993l ) proposed that three effects could combine non- 
linearly to inhibit the detection of CN and C2 in Comet Yanaka. First, cosmic ray exposure 
could "glue together" or "carbonize" the comet (deplete H, N, and O relative to C), with two 
effects: trapping carbon-bearing gas molecules between dust grains, and producing larger 
dust grains more resistant to heating and fragmentation. Lastly, since Comet Yanaka was 
observed over a smaller spatial dimension (660 x 16,470 km) than the other comets in Fink 
(1992; 3420 x 84,200 km and 2400 x 145,600 km), there may be less time for molecule 
production to occur within the spectrograph slit aperture. 

Our observation of 96P covers a spatial area of about 820 x 7 0,570 km. Our sl it length 
is comparable to the larger spatial observations in iFink I (119921 ). I Cochran I (119851 ) reviews 
Haser-model decay scale lengths for the parent molecules of CN (1.2-2.2x10^ km), C2 (1-7- 
3.5x10^ km), and C3 (1.0-2.5x10'^ km), which are all shorter than our slit length. If cosmic 
ray gluing increased the parent decay scale length of C3 by an order of magnitude, we 
still expect to see C.s production. Brightness and polarization measurements of dust in the 
coma of 96P (IGrynko. Jockers. fc Schwenn 1 12004| ) show grain sizes consistent with in-situ 
measurements of Halley dust, also pointing to limited cosmic ray gluing. 



Greenberg. Singh, fc de Almeida I (119931 ) admit that other dynamically new comets do 



not display extreme carbon depletions. Likewise, other short-period comets {P < 7 years), 
presumably with surfaces similarly processed b y solar radiation as 96P, display measurable 
ca rbon production: 2P /Encke and 6 P/d'Arrest (lA'Hearn. Millis. &: Birch Ill979l ) , 9P/Tempel 



1 ( Lara et al. l2006h . 1 9P/Borrellv (IHamane et al 



1993h. 46P/Wirtanen (| Jockers. Credner. fc Bonev 



20021 ). 26P/Grigg-Skjellerup (IJockers et al 



19981 ). and 67P/Churyumov-Gerasimenko 



Schulz. Stiiwe. fc Boehnhardt II2004J ). In old "dusty" comets, or comets moving outbound 
from perihelion, the CN and C2 emission bands are the last to be observable as global gas 
production turns off (i.e. A'Hearn et al. 1995, Fink & Hicks 1996), such that CN is used as 
a gas-to-dust indicator (i.e. Storrs et al. 1992). However, CN is not detected in the Kast 
spectra of 96P. 

Other dormant comets may perhaps show mostly ices, and n ot carbon-b e aring gas, when 



they outburst. However, CN was det ected in 9 5 P/Ch iron by iBus et al. 
4P/Faye was observe d in outburst by Grothues (Il996 



C3 production rates (iGil-Hutton fc Licandro I Il994l ). 



(Il99lh . Comet 
and has measurable CN, C2, and 



Chamberlin et al. I (119961 ) detected 



no CN in three candidate comet-asteroid transition objects, although their signal-to-noise 
deteriorates in the region of the CN band. Unfortunately, none of these observations were 
programmed to look for the NH 3360 A band. 



Having exhausted several explanations for the odd spectrum of 96P, we can compare it 
to other comets with extreme production-rate ratios. The values in square brackets in Table 



- 6 - 



m refer to logarithmic ratios between the production rates of variou s species. Wher e there 



are discrepancies between the [X/OH] and the [X/CN] ratios for the lA'Hearn et ah I (119951 ) 



comets, the [X/CN] ratios are more rehable, and are the basis for the hsted [X/NH] ratios. 
In 96P, we confirm the S07 report that the carbon-chain molecules are depleted compared to 
NH. Even compared to other "C-chain-depleted" comets 96P is unusual because we detect 
C2 emission but not CN. A distinctive property of 96P is the very high [C2/CN] ratio, for 
which we only obtain a lower limit. Rather than being carbon-chain depleted relative to CN, 
96P is the comet with the least amount of CN to also display a C2 band, and by far the comet 
most depleted in CN relative to NH, C2, and C3 than either "typical" or "C-chain-depleted" 
comets. 



5. CONCLUSIONS 

Comet 96P/Machholz has extremely depleted C2 and C3 relative to NH and NH2, with 
CN more extremely depleted. The dimensions of our spectrograph slit were sufficient to see 
any possible decay of the parents of these carbon molecules. The extreme carbon depletion of 
96P is unlikely to be due to cosmic ray "gluing" in the comet's past, and unlikely to be due to 
surface processing by the Sun over repeated short orbits (since this effect is not seen in many 
other short-period comets). Though 96P was observed to have an outburst shortly after 
its discovery, other confirmed outbursting comets show carbon emission features. Thus, 
it appears that 96P/Machholz belongs to a small class of comets with genuinely unusual 
molecule production rates. 

This research is supported by NASA grant NNG05GG59G through the Planetary As- 
tronomy program. 
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Table 1. Comet 96P/Machholz: Observations and Orbit Parameters 



UT 


Th (AU) 


A (AU) 


Airmass Exp (s) 


2007 Apr 27 11:38 


0.754 


0.747 


2.5-1.87 2 X 1200 


q = 0.1241>i;7 


e = 0.9588 




i = 60.1813° 
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Table 2. Haser Model Scale Lengths and Gas Production Rates Q for Comet 

96P/Macliholz 







Spectrum 1 






Spectrum 2 




Mean 


Molecule 


Parent 


Daughter 


Q 


Parent 


Daughter 


Q 


Q 




(10^ km) 


(10^ km) (10^5 


mol. s~^) 


(10^ km) 


(10^ km) (10^5 


mol. s~^) 


(10^5 mol. s-i) 


NH 


70 


200 


10.54 


40 


120 


6.18 


8.36 ± 2.18 


NH2 mean 


3.4 


42 


26.22 


4.2 


25 


33.54 


29.88 ± 3.66 


-(0,10,0) 


3.0 


50 


17.07 


6.0 


25 


23.25 


20.16 ± 3.09 


-(0,11,0) 


3.0 


40 


20.3 


4.0 


25 


21.65 


20.98 ± 0.68 


-(0,12,0) 


5.0 


35 


13.03 


5.0 


>300 


29.15 


21.09 ± 8.06 


-(0,13,0) 


2.0 


160 


47.14 


3.0 


150 


64.47 


55.81 ± 8.67 


-(0,14,0) 


4.0 


>300 


33.56 


3.0 


210 


29.19 


31.38 ± 2.18 


CN 


7.5 


240 


NA 


7.5 


240 


<0.0075 


<0.0075 


C2 


7.0 


20 


0.039 


7.0 


30 


0.051 


0.0452 ± 0.0061 


C3 


1.3 


57 


<0.020 


1.3 


57 


NA 


<0.020 
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Table 3. C/1998 Yl: Observations (Fink 1992) and Orbit Parameters 



UT 


rn (AU) A (AU) 


Airmass Exp (s) 


1989 Jan 15 13:26 


0.932 0.367 


1.84 3 X 300 




e = 1.0000 


i = 71.0077° 
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Table 4. Comets With Extreme Production Rate Ratios 



Comet 


Type 


[CN/OH] 


[C2/OH] 


[NH/OH] 


[C2/CN] 


[CN/NH] 


[C2/NH] 


[C3/NH] 


Ref. 


"Typical" mean 


T 


-2.50 


-2.44 


-2.37 


+0.06 


-0.13 


-0.07 


-1.22 




"C-Chain Dep." mean 


D 


-2.69 


-3.30 


-2.48 


-0.61 


-0.21 


-0.82 


-1.70 




P/Howell 


JF,T 


-2.47 


-2.55 


-1.14 


-0.09 


-1.33 


-1.42 


-2.68 




Bowell (1980b) 


DN,D 


-3.33 


-3.43 


-2.06 


-0.39 


-0.58 


-0.97 


-2.26 




P/IRAS 


HF,D 


-2.57 


-3.48 


-2.03 


-0.93 


-0.45 


-1.38 


-2.28 




P /Wolf-Harrington 


JF,D 


-2.53 


-3.75 


-2.78 


-1.20 


-1-0.30 


-0.90 


-1.75 




Yanaka (1988r) 


DN,D 


<-3.97 


<-4.34 


-3.14 




<-0.83 


<-1.20 




2 


96P/Machholz 


HF,E 


-4.93 


-3.63 


-1.86 


-1-1.30 


-3.07 


-1.77 


-3.17 


3 


96P (this work) 


HF,E 








>-|-0.78 


<-3.04 


-2.26 


<-2.61 





Note. — Comet families: JF=Jupiter, HF=Halley, DN=Dynamically New. "Carbon-chain depletion" relative to CN: 
T=Typical, D=Depleted, E=Enhanced. Very unusual ratios are highlighted in bold-face. 

References. — (1) A'Hearn et al. 1995; (2) Fink 1992; (3) Schleicher 2007 
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Fig. 1. — Atmospheric extinction in magnitudes per airmass versus wavelength at Lick 
Observatory on 2007 Apr 27 UT. The sohd hne is the measured extinction. The long-dashed 
line is the ozone component. The short-dashed line is the Rayleigh scattering component. 
The dot-dashed line is the aerosol component. The crosses show extinction values tabulated 
previously for Lick Observatory. 
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Fig. 2.- 



erg s ^ cm ^ 



Average integrated spectrum of 96P. (a) Flux per unit wavelength, scaled to 10~^^ 
Note the strong NH line at 3360 A. (b) The same averaged spectrum, 
scaled to 10"^"^ erg s"^ cm~^ A~\ to clearly show the NH2 series from 4500-5800 A and the 
faint C2 band at 5165 A. The CN band and the C3 band at 4050 A are labeled in parentheses 
to denote upper limits. The questi on marks refer to lin es th at are unidentified in the hig h- 
resolution comet spectral atlases of iBrown et al. ] Jl996h and I Cochran &: Cochran I (120021 ). 



